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Abstract . In this paper we have considered the problem of the 
structure of a strong collision free wave in a plasma perpendi- 
cular to its magnetic field. It is shown that conditions on 
the two sides of the wave are the same and are separated by a 
solitary pulse provided the kinetic energy of the motion at one 
end of the wave (x = say) is less than the sum of magnetic 
and thermal energy of the plasma. If the kinetic energy is more 
than the sum of magnetic and thermal energy, no wave exists aruj 
the constant solution at x = -° is the onl y solution of the 
equations describing the wave form. 
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I. Introduction. 

The subject of the structure of a strong collision free wave, in a 
plasma perpendicular to the m a g netic field has been treated by Adlam and Alan\ 
Burgers^ and Haiij Lust and Schluter^. Hain, Lust and Schluter in particular 
have obtained solutions of the wave structure by solving the non-linear 
equations on a computer. 

In this paper we have investigated mathematical properties of the equa- 
tions which help us to guess the correct solution out of the many solutions 
obtained on a computer. As an illustrative example we shall consider a 
plasma with isotopic pressure. The generalization to the anisotropic case 
is qpite straight forward. Burgers has made an extensive study of the wave 
structure when the plasma pressure is anisotropic and he arrives at the same 
general result that there is no change of state on the two sides of the wave 

1. e., that no shock wave exists in a collisionless plasma perpendicular to 
the magnetic field. 

The result of our paper is that a solitary wave only exists if the 
kinetic energy of motion, which produces the wave, is less than the sum of 
magnetic and thermal energy of the plasma; that no stationary wave exists if 
the kinetic energy of motion is more than the sum of magnetic and thermal 
energy of the plasma. In the latter case the constant solution at one end of the 
wave is the only solution of the differential equation (section 4). 

2. Equations Describing the Wave Form. 

Consider a stream of plasma at x = -“ with a constant velocity u & along 
x-axis and a constant magnetic field B^ along z-axis . The electrons and 
ions of the plasma will drift in the y-direction because of the magnetic 


field 
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Our object is to obtain conditions at x = + 00 and the form of the wave 
or waves which start from the given solution at x = -®. As a simplifica- 
tion we assume that the electrons and ions satisfy Maxwellian distribution 
function with variable number density, velocity components and the tempera- 
ture. Therefore, if f e , fi denote the electron and ion distribution 
functions 





M 



) 



In cartesian coordinates, assuming steady state, the Boltzmann equa- 


tion can be written as 
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and similar equations for the ions . Let us call the similar ion 
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ixvoll's equations yield 
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To further simplify the problem, we will make the assumption that 
the space charge can be neglected. The characteristic length 

A 

of space charge oscillations is of the order of the Debye lengtlr 

/ kT 

{? = [ ) while, as we .shall see in this paper, that of magneto- 

\4Kney . 

hydrodynamic oscillations is of the order of J itnUe % * 

Thus assuming that the electron thermal velocity is well below the 
velocity of light, the space charge oscillations are of much lower 
scale than the oscillations we are interested in. With the as- 
sumption k? 0 and —"*1 (to a good approximation 


we then have 


= u^. — a > say), equation (6), (6)' and 


(8), (8)' yield, respectively 

J-[« n u^^kT c ^kT e ] = -J. 
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Equation ( 5 ) now gives 


>UA - **>» 


(15) 


where Limit n(x) = n-, and Limit u(x) = ui . 
Equation (13) can ho integrated and yields 


R 1 1 
M'nu? * ’'kT^ + ~~ — _ M-»,u*+27>kli+ t 
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where Limit B (x) = B . Eliminating T.^. T from equations (l 6 ) 

X-^-00 1 le 

and (14) we obtain 
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■M t 1 


=0 
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Dividing by 


J <frr My^uj" 


and using the transformations 



I _ u 
J 8 -k 


kJlaJt _ n r' _ 

Mu, 1 


we finally obtain the equations 



3i e, 
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and 


(19) 


( 20 ) 



where 
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l + <2 Oj •+■ 



( 21 ) 


Equations (19), (20) satisfy the boundary conditions H = jb k 5 U 1 I <s^ J 


.*0 


3* A Mathematical Discussion of Equations (19), (20) 

The pair (U, H) = (l, |:!.j ) is a constant solution of these 
equations. We shall study other constant solutions represented 
by (19), (20), but we first examine the stability of the solutions 
of (19), (20) in the neighborhood of (l, and determine the 
nature of this critical point . 

Let U = 1+U^, H = We then obtain from (19) 

H, = U, F ( H, d Uj) 

where 

( s h + H.)0 + vJ 


( 22 ) 


( 23 ) 


From (22) and (23) we obtain 


- 7 - 


where 


L „ FM = (25) 

4* » F(°><0 ^ M * ^ (•'*) 

dH* dUi (26) 

= _ Q^i-p 1 <fo 

etc. s- /a, 3 p. 

Equations (20), (24) yield, 

g. *■ y,(u,) t k o, ♦ -• (27) 

where ^ ^ 7^ are analytic in near =0, ^(U^) = 0(U^ 2 ) and 


1/ _ , P _ i- ~ fh 

/S. - H r- =* 7 

>-» 1-601 

To determine the nature of the point U = 1, we linearize equa- 
tions (23), ( 27 )* The linearized form of equation ( 27 ) is 


( 28 ) 


1 + K», 

In all the discussion that follows we shall assume the flow 
velocity to "be much higher than the thermal velocity of particles 
constituting the plasma, so that 60-^ 1 or 1 - 60^^ 0. The sign 



( 29 ) 
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of k^ is therefore determined by the expression 1-2 ft 2 ~ ^ 9 1 ‘ 

In case I when i - 2 A 2 - *x> 0 , k 1 is positive and therefore 
U = 1 is a center in both (27) and (29) according to a known result 

of Lyapounov, see Malkin , pp 123-4, (as no odd powers of dU^/dy occur). 

In case II when 1 - 2 - 60 ^ 0,1c 1 is negative; the point U = 1 

is a saddle point in both equations ( 27 ) and ( 29 )* 

Let us now consider all the constant solutions represented 
by the equations ( 19 ), (20). According to equation (20) constant 
solutions satisfy the equation HU = j$^. This condition, together 
with equation ( 19 ) then gives the equation 

/CO) H U* 1 - uVt‘ 4,i *' + l 3 ^) + / Sj ^ '+3-04=0 ( 50 ) 

f(U) may also be written as 

fOi) = (V<) [<J 3 - S) - 4 9 ‘] ( 51 ) 

- (U" 1 ) fXu) 

From Descarte's rule of signs h(U) cannot have more than one 
positive root. Also it is easy to see that it has a positive root. 

Let this positive root be denoted by Uo* Since f(0)^0, Uo will be 
less than or greater than j. when f'(l)^> 0, f'(l)< 0 respectively. We 

therefore again have two cases: 
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(i) f(i) = 1 - 2 y&j 2 - 6e 1 > o 

and 

(II) l - 2 pf - 6e 1 <0. 

As we require that 60 1, Case I may occur. Figure l shows 

a rouch sketch of f(U) versus U for cases I and II above. 

To find the nature of the other critical point U = U 0 , we 
write equation (20) as 

Jy /rfyfr vftu)+ £Cj )| 

[<%) l u p'Lu) J + ° <32) 

where 

r .. 2. 0 1 1 ^ 

H= = [ A ~ U ~ v z] (55) 

g'(U) has a zero at IP = It is interesting to note that 

this si n gularity of the differential equation can be quite close to 1. 

For example, for 0- L = .1, v* .880117. Higher values of 0 1 

can bring this singularity still closer to U = 1. 

When equation (32) is linearized about the zero Uo of f(U), 
we obtain the equation 

u, +Ua) 0 , 


=. o 


( 34 ) 
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where 0 = U 0 +U, and U L (U) - 1 + ? r 1 . The nature of 

U s (6’.-U*) 

the point Uo = 1 has already been discussed. Now let us assume 
in addition to 69 ^ 4, that Uo 4 ^ 60^. We wish to establish 
for equation (3a) that in ease 1 when the point U ■ 1 is a oenter 
the other equilibrium point Uo 41 is a saddle point and conversely, 
in case II when U = 1 is a saddle point, the other equilibrium point 
Uo ^1 is a center. By arguments similar to those given above for 
equation (27) and (29), this is equivalent to showing that if 
h(Uo) = 0, h given in (31), then Uo^l implies L(Uo) ^ 0 and Uo 1 
implies L(Uo) ^ 0. 

It can be shown that 


U?Cu)- /*, 


- fiv) 

ulp, + v? CO)] 


(35) 


Now, consider case I. f (U) has a zero at U = Uo (4 l) and U = 1. 

In the U,U plane construct a closed loop C enclosing Uo and 1 and 

• # 

lying so close to the U-axis that the sign of U is determined by 
the last term on the left in (32). (see Figure 2). Then using 
(35) and upon considering the changes, in the sign of f(U), it 
follows that the index of C with respect to (32) is 0^. As the 
index of a center (U = l) is 1, this implies that the index of Uo 
is -1. If L(Uo) f 0, i.e., if Uo is an elementary critical point, 
then Uo of necessity must be a saddle point and L(Uo) £*.0 will hold. 

A similar arrangement for Uo^l in case II shows that the index of Uo 

is 1 and therefore L(U 0 ) > 0 holds (provided of course L(U 0 ) f Cli 
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Thus the problem is reduced to showing that the two poly- 
nomials h(U) = U 3 - bU 2 - b'J - 20, and M(U) = U 6 - 6Q ± U 2 *^ - 2b 
have no common positive zeros in either case I or case II. Here 
we have used the definition b = 20 ^ 2 in deriving the form of 

M(U) which is the numerator of UL(U). 

This has not been shown in general, although for particular 
choices of ) 0-^, it is relatively straightforward. The critical 
case, however, is when b + 0^ = ^ or 2 A 2 ^. 60 ^ = 1, for then the 
two roots Uo and 1 become equal and, of course, in this case M(l) = 0. 

Let us suppose that b+0, is close to X i.e., Uo is close to 1 
so that (Uo - l) 2 , (Uo - l) 3 are negligible in comparison with Uo - 1. 

Now let R = U 2 and expand h and M about U = 1 and R = 1. This gives 

&(v) = i - zck+oi) + 3O-^)C^~0 +(s-k)(u-i?+(y-!f (36) 


and 


M(K) = 





The zeros Uo, Ro (=U 2 - l) of 

Uo =r | + 

R. — » + 


and M may be approximated by 

SO- 

3 6-2&t) 


( 32 ) 
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Now suppose we are in case I, that is b + 0^ ~ . Then as 

we have Uo - 1<R 0 - 1 = U 2 - l<Oj where M(U) =0, U >0. This 

implies 

i-0„ > ',-v)Ci+ OjrO-v) 

Therefore 

U e < U 4 I 

1 — 

In case II, b+ 0 .^— , it follows similarly from Uo - l^U 2 - 1, 

that 1<U ^U 0 . Thus the positive roots of M and R, which coincide 
for b -i- 0 1 = do separate as required, for b -(- 0 ^ near Although 
a proof for the general case is lacking, the following calculations 
may be of some interest. 

Taking Q ± = .1, b = .3, we find (60 1 ) 1// ‘^ = .880117, U 0 — .88978 
and M( T Jo) = - .1788 <0. With 0^^ = .1, b = . 5 , we find U 0 as 1.1116 
and M(Uo) ^ -5^567 ^ 0. Thus we observe that at b = . 3 , Uo is very 
close to U* = (60^)^"' ^ and so b could not be decreased further and 
still have Uo 4 ^ 60 

Thus we have seen that in case II, that is Uo ^1, if M(Uo) > 0 , 
then Uo is a center. In this case the phase plane portrait of the 
solution is shown in figure 3- In the other case, in which Uo is a 



saddle point and 1 a center, the phase plane picture is the same as 
figure 3 with the role of 1, Uo reversed. As is clear from our 
analysis we have tried to determine the structure of the solutions 
of equations ( 19 ), (20) from the properties of the differential 
equations (27) or (32). In the next section, we shall verify our 
predictions by' obtaining some solutions on a computer. 

4. Computer Calculations 

Some of the results given in the previous section were veri- 
fied by computing the solution of equations (19) , (20) on an IBM 
7094. The way to do this is to choose values of U, H (different 
from 1, lying on the eigen solution by using the linear equation 
(29). All calculations were started close to U = 1, H = Since 

the solution depends critically on the position of the point of 
start in the phase plane we ran several cases using values slightly 
different from those obtained from the linear equation (29)* For 
example, for ^ = .81, 0 1 = .05, we obtained the solution using (a)U=1.0001, 

= .00001 and (b)U=1.0001, = .001 while from the linear equation 

ay ay 

we obtain a -Blue of ^ which lies between .00001, .001. Figures 4a and 4b 

dy 

give the solution for the casea (a), (b) above. In the case (b) 

dU 

the solution becomes unstable because we happen to use a value of ~ 
which lies outside the oval of figure 3- In the case (a), we obtain 
a solution which almost returns back to its initial value and re- 
peats itself. (The repetitions are evidently the cause of the computa- 
tional error, which if decreased will show that the solution is a 
solitary pulse. Since only two branches of the curve approach 
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U = 1, U = 0 the solitary pulse is the required solution of the 
wave form. ) 


We now turn to a calculation of case I considered in section 3- 
We used ^ = .1, 0^ = .1. Taking a point away from U = 1, “ = 0, 

we obtained a solution consisting of oscillations having U = 1 


= 0 for its center in accord with the fact that U = 1, -t— = 0 

* dy 

is a center. We see therefore that in case (i) even the solitary wave 


does not exist. 


We have not made any computations about the point U = Uo as this 
point is of no interest to us. 

This corresponds to a solution lying just within the oval pictured 
in figure 3- A more accurate choice of initial conditions would show 
that the boundary of this oval corresponds to a solitary pulse, that 
is, the solution forming the boundary of the oval leaves U = 1, U = 0 
at y = - co, and returns to this point at y = -t- oo . 
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Figure Captions 


Fig. 1 . f ("y) vs U for Cases I, II of text. 

Fig. 2. The closed curve C enclosing the points U = U and U = 1 
and lying close to the U axis. 9 

Fig. 3 . Approximate phase plane picture of the solution for the 
case II when \J = 1 is a saddle point and = Uo is a 
center. 


Fig. 4 a. Computed solution of equations (19), ( 20 ) for 

6-l = .05, using perturbed values (J = 1.0001, dU _ 

<3y " 


= .8l, 
. 00001 . 


Fig. 4 b. Computed solution of equations (19), ( 20 ) for 

e x = .05, using perturbed values TJ = 1.0001, dU 

ay 


. 001 . 
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Fig. 4b 




